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Chirped Dissipative Solitons 

Vladimir L. Kalashniko^fl 
Institut fiir Photonik, TU Wien, Gusshausstrasse 21/381, A-IO4O Vienna, Austria 

The analytical theory of chirped dissipative soliton solutions of nonlinear complex 
Ginzburg-Landau equation is exposed. Obtained approximate solutions are easily 
traceable within an extremely broad range of the equation parameters and allow 
a clear physical interpretation as a representation of the strongly chirped pulses in 
mode-locked both solid-state and fiber oscillators. Scaling properties of such pulses 
demonstrate a feasibility of sub-mJ pulse generation in the continuous-wave mode- 
locking regime directly from an oscillator operating at the MHz repetition rate. 

PACS numbers: 42.65.Tg, 42.65.Re, 42.55.Wd, 42.65. Sf 
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I. INTRODUCTION 



Q I The nonlinear complex Ginzburg-Landau equation (CGLE) has such a wide horizon of ap- 



plication that the concept of the world of the Ginzburg-Landau equation has become broadly 
established [l]. The CGLE is used in quantum optics, modeling of Bose-Einstein condensa- 
tion, condensate-matter physics, study of non-equilibrium phenomena, and nonlinear dynam- 
_0 I ics, quantum mechanics of self-organizing dissipative systems, and quantum field theory. In 
^ ' particular, this equation being a generalized form of the master mode-locking equation 0, 0] 
r-| ■ provides an adequate description of pulse localized inside a laser oscillator or propagating along 
Ph. a nonlinear fiber. 

The exact analytical solutions of CGLE are known only for a few of cases, when they 
CN ' represent the dissipative solitons and some algebraic relations on the parameters of equation 
are imposed Jsj. More general solutions can be revealed on basis of the algebraic nonperturbative 
techniques j^, which, nevertheless, are not developed sufficiently still. 
^ I The perturbative methods allow obtaining the dissipative soliton solutions in the vicinity 

of the Schrodinger solitonic sector |5| of the CGLE [6j. Another approximate method utilizes 
^ ■ the reduction of infinite-dimensional space of CGLE to finite-dimensional one on basis of the 
O ! method of moments 01 ■ 

In this article I present the direct approximate integration technique for the generalized 
CGLE [sl, 0] . This technique is suitable to a physically important sector, which is represented 
^ ■ by the chirped solitary pulse solutions, or the chirped dissipative solitons (CDSs) of the CGLE. 
The method under consideration does not impose any functional constraints on the solution 
and equation parameters. Additionally, the dimensionality of the obtained CDS is reduced 
in comparison with the parametric space of the CGLE. As a result, the CDS characteristics 
become easily traceable. 



As was found [lOj, the CDS is energy scalable due to strong pulse chirp, which results in 
the substantial soliton stretching. Such a stretching reduces the CDS peak power and, thereby, 
leads to the soliton stabilization. Simultaneously, a strong chirp provides a sufficiently broad 
spectrum allowing the pulse compression by the factor 10-=-100. As a result, one can obtain the 
femtosecond pulses with over-/iJ energy at the fundamental MHz repetition rate of an oscillator. 
This allows the table-top high-intensity experiments with the level of intensity of the order of 
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10^^ W/cm^ 28|. Thus, an analysis of the CSP solutions of nonlinear CGLE is of interest not 



only from the theoretical but also from the practical points of view. 

II. APPROXIMATE INTEGRATION OF THE CGLE 

Let the 1+1- dimensional generalized CGLE be written down in the following form j^, 0] 

d (9^ 

—a iz, t) = -aa (z, t) + a—ra (z, t) + \\a (z,t)\ la (z, t) + 

OZ Ot^ /IN 

-i {'y + x\a'{z,t)f) \a {z,t)f a {z,t) , 

where the field envelope a{z, t) depends on the propagation (longitudinal) coordinate z and 
on the local time t ( "transverse" coordinate; for an interpretation of ([1]) in the terms of Bose- 



Einstein condensate see, for instance, |12l . Il3|). The real part of the right-hand side of ([T]) 



(nonlinear Ginzburg-Landau equation) describes the dissipative effects: i) net-loss with the 
coefficient a, which is the difference between unsaturated net-loss and saturated gain; ii) spectral 
dissipation with the squared inverse bandwidth coefficient a, and iii) nonlinear gain, which is 
some function E of the power |ap. The imaginary part of the right-hand side of ([T]) (nonlinear 
Schrodinger equation) describes the nondissipative effects: i) group-delay dispersion (GDD) 
with the coefficient /3, and ii) self-phase modulation (SPM) with the coefficient 7. The SPM 
can be saturable or enhanceable with the corresponding coefficient x- 

The nonlinear gain (or self-amplitude modulation (SAM)) in Eq. ([1]) describes the depen- 
dence of dissipation on the field. To provide a soliton formation, E has to grow with power at 
least for a small P{z,t) = |a(2;,t)p. In an oscillator, the SAM can be provided by the different 
mechanisms: self-focusing inside an active medium, nonlinear polarization dynamics, bleaching 
of an absorber (e.g., a semiconductor quantum- well structure), etc. (see [l, [13] )■ Below two 
main functional forms of S will be considered: i) cubic-quintic SAM term E = kP — (kP^, and 
ii) perfectly saturable SAM fiKP/{l + KP). Here k is the effective inverse loss saturation power, 
( is the parameter of the SAM saturation >0 to provide the stability against collapse), and 
/i is the SAM depth. The first type of SAM describes a Kerr-lens mode-locked oscillator, and 
the second type corresponds to an oscillator mode-locked by a semiconductor saturable ab- 
sorber. It should be noted, that the functional form of E is not decisive for the method under 
consideration. 

The main assumption underlying our method is that the dissipative soliton has a large chirp. 
Mainly, this assumption is valid in the positive dispersion range (/3 >0) of Eq. ([1]), but the is the 
CDS also in the negative, or anomalous, dispersion range {(3 <0, see [l3|)- As a result of large 
chirp, the CDS is temporally stretched that admits the adiabatic approximation (T ^ a/5, T is 
the soliton width), and the fast variation of soliton phase admits the stationary phase method 
in the Fourier domain (igI . 17| . Formally, the approximations under consideration mean that 



/3 ^ a and 7 ^ k (or 7 ^ fxn). That is the nondissipative effects prevail over the dissipative 
ones. 



A. Cubic-quintic CGLE 

Let the SAM term in Eq. ([T]) be written in the form 

E[|a(^,t)|'] =KPiz,t)-(KP{z,tY, (2) 



2 



Chirped Dissipative Solitons 



and the desired travehng wave solution of ([T]) be 



a {z, t) = a/ P (t) exp (t) — iqz] . (3) 

Here (pit) is the phase and q is the sohton wavenumber. Substitution of (12]) and (E]) in ([1]) with 
subsequent apphcation of the adiabatic approximation results in jo], [isl, [l9| 



^,'ndp dn\ . , . 



(4) 



where Q = d(l){t)/dt is the deviation of instant frequency from the carrier one. The latter is 
assumed to coincide with the minimum of spectral dissipation. The solution of (j4j) is 

P = ( + 4gx - ^xm^ - 7)/ 2x, (5) 
which exists in the limit of x ~^0- 

P={q-Pn'')/^. (6) 

The condition of power positivity requires an existence of maximum frequency deviation 
= q/P. 

Second equation in (jl]) leads to 



dn 
It 



a + aVL^ - 



2 _ >t(A-7) 



4^2 



(2x + C7 - C^) 



P[Axm^-{A -^)A] 



(7) 



In the limit of x ~^0) has |17|, [20 



2 37 / . c . //. c 



7P(0) = ,3A^ = ^(l--±y^l-- 



4(tC 



7 



/3S^ = ^(l + c)--7P(0), 

where c = a'yj fin and the equation for the peak power P(0) (or, equally, for A) results from the 
regularity condition dVL/dt < ±00. Such a condition for Eq. ([7]) leads to [q = a( / k, rj = / x)'- 



7 



16C/3 + 1) 



2(3+^ + ^ 



c I 3r? _|_ 
2 2 



2f-16, (1 + 1)) 



- 12 - 3c - 9r7 - — 

7] 



(9) 
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In the cubic nonlinear limit of Eqs. (ITP]) {x — j-O, ( — j-O), which admits an exact CDS solution 
(2l| . one has 



^ c 

_ 37 y _ 37 

^ /t(l + c)' fi;A(l + c)' 

One may see that the maximum parametric dimension of (|T|2|) is 4, but the parametric 
dimension of the approximate solutions is reduced (the maximum dimension corresponds to 
(l5|7Il and equals to 3, i.e. rj and c). 

Next assumption allows further simplification. Since the phase (j){t) is a rapidly varying 
function due to a large chirp, one may apply the method of stationary phase to the Fourier 
image e{u) of a{t) [l6|. Then, the spectral profile corresponding to ( 1S|7|) can be written as 



„ vr M - 1) (M - 1) cr7 + 4 (2u^ - A^)) H (A^ - 

p[u) = \e{u) - 



cA {{A -l){c{<^ + r] + r]^ + uj^)+r] {A^ - co^)) - 2 (r/ + 1) (A2 - 00^)) ' 

(11) 

where the normalizations: t' = t{K/() ^ KjaC,, A'^ = A'^aC^/K^ VL'"^ = Q'^aC/n, P' = (P, and 
E' = E (k/'j) -sJ K(/a are used (primes are omitted). A = a/1 + 4 (A^ — Vfi)lrjc and H is the 
Heaviside's function. 

As it has been demonstrated in [9!] , the truncated at ± A spectra (11 II) have convex, concave, 
and concave-convex vertexes (Fig. [1]). When x ^0, only convex truncated spectra remain |17| : 

, , 67r7 H (A2 - u;2) 




FIG. 1: Spectra corresponding to the — branch (i.e. — sign in Eq. Q) of CDS for c =1, ? =0.01 and 
different values of tj. t] >0 corresponds to an enhancement of the SPM with the power growth and 
\r]\ — 00 corresponds to an absence of the quintic SPM, i.e. \x\ -^0- 

At last, for the cubic nonlinear limit of Eq. ([1]) {x — j-O, C — t-O) one may obtain 
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That is the truncated flat-top spectrum, which approximates the spectrum of exact solution of 
the cubic CGLE ([T]). However, the exact spectrum is expressed through the beta-function that 
troubles its analytical study. 

An appearance of the concave spectrum f|T2l) from the flat-top one f|T3|) can be illustrated 
by means of the perturbation analysis in the spectral domain. The approximate method under 
consideration allows easily performing such an analysis. In the case of the cubic CGLE with 
the real quintic term as a perturbation, the equation for the Fourier image of the soliton 
perturbation / {u) can be written as [3| 

oo oo 

[k (u) - q] f iu) + ^ J du'U (u - u') f (u') + i- y du'V (u - u') f (c^') = S (u) , (14) 

— oo — oo 

where k (u) = —/Su"^ — i{a + au'^) is the complex wavenumber of linear wave, the kernels 
U and V are the Fourier images of and a(t)^, respectively (a(t) corresponds to the 

time- dependent part of solution ffTOl) and co-propagation of a stable CDS with a sustained 
perturbation is assumed). The assumption, that the perturbation is in-phase with the CDS 
allows the replacement ® /* -H- ?7 ® / in ffl^. It should be noted, that such an assumption is 
only conjecture for the CDS (compare with [2^). The Fourier image S{uj) of the perturbation 
source can be obtained from the perturbation term —K(P(tya{t) by means of the method of 
stationary phase (see above): 



:i + c) 



2k {{A'' - u^) {1 + c)) 



(15) 



For the kernel U, one has 



U{uj- u') 



97r7/3 (7 + ik) 



oj') csch 



2kA (1 + c) 



(16) 



Resulting equation is the Fredholm equation of second kind and can be solved by the Neu- 
mann series method so that the iterative solution is 



fn (w) 



S{u) 3 
k {u) — q 271 {k (cu) — q) 



du'U {Uj-u')fn-l{uj') 



(17) 



where /„ (w) is the n-th iteration and /o (w) = S {oj)/ {k {u) — q). 

The successive terms of f[T7|) are shown in Fig. |5J One can see that the tabletop spectrum 
becomes convex and some oscillating substructure develops. The last phenomenon explains 
the results of 23|, where the cause of both regular and chaotic pulsations of CDS has been 



attributed to an excitation of the solitonic internal modes revealed by oscillating structure in 
Fig. El 

It is convenient to represent the CDS existence regions in the form of the so-called master 
diagrams 17|. The coordinates of such diagram are the normalized energy and the c-parameter. 
The energy E = (u) can be easily obtained from Eqs. (Illfl2fl3p and connected with 

the (T-parameter through a ~ 6{E/E* — 1), where 6 = da / dE\^^^^ and E* is the energy of 
continuous- wave solution of Eq. ([1]) [l7l |. 

The master diagram for the CDS is shown in Fig. [3] for the case of vanishing quintic SPM 
(x =0)- The CDS is two-parametric in this case. The solid curve shows the border of the 
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i,u ' r - . ' r — ' 1 ' 1 ' 1 ' — ^ ' ^ ^ 1 ■ 

-0,015 -0,010 -0,005 0,000 0,005 0,010 0,015 

a, fs"' 



FIG. 2: Absolute values of the successive terms in the Neumann series ()17p for the real quintic nonlinear 
perturbation (tildes are omitted), a =16 fs^, f3 =250 fs^, k = 0.047, a = 2.5 x W~'^, and ( = 0.2j. 

CDS existence q =0. Below this border, the vacuum of Eq. ([T]) is unstable. The gray curve 
divides the existence regions for the + and — branches (the corresponding signs in Eq. dH])). 
The branches merge along this curve. Crosses (circles) represent the curve along which there 
exists the +(-) branch for some fixed value of ? (isogain curve). 

o 
o 




1 ' I'o ' 100 



E 

FIG. 3: The master diagram for the CDS of Eqs. ()ll2p . x =0- There exists no CDS below the black 
solid line (here a <0, i.e. the vacuum of ([1]) is unstable). Black solid curve corresponds to =0. Gray 
curve divides the regions, where the + and — branches of ([9]) exist. Crosses (circles) correspond to 
the +(— ) branch for ^=0.01. 

The master diagram reveals four significant differences between the + and — branches of 
CDS. The first one is that the — branch has lower energy than the + branch for a fixed 1/c. 
The second difference is that the + branch of an isogain has a finite bound for E ^ oo. In this 
sense, the + branch is energy scalable, that is its energy growth does not require a substantial 
change of c. The — branch is not energy scalable, that is the energy growth along this branch 
needs a substantial decrease of c (e.g., owing to a dispersion growth). The third difference is 
that the + branch verges on ^=0 within a whole range of E. The fourth difference is that the 
— branch has a Schrodinger limit X (Eq. (11011131) ). 

For X 7^0, the CDS is three-parametric. Contribution of the positive quintic SPM (i.e., x >0 
and the SPM enhances with power) narrows the CDS existence region (Fig. S]). This means 
that smaller c is required to provide the CDS existence for a fixed E. That is, since the positive 
quintic SPM means an enhancement of the SPM with power, a SPM enhancement has to be 
compensated, for instance, by a dispersion increase c oc 1/(3. One can see from Fig. HI that the 
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+ branch region narrows substantially with a growth of positive x within a whole range of E. 
However, it is important to note, that the range of c, where the CDS with a fixed <j exists, is 
defined by the difference between i) the point of intersection of the isogain with the boundary 
between the + and — branches and ii) the isogain asymptotic iov E ^ oo. As a result, the 
range of c, where some isogain exists, broadens with x >0 (Fig- HI)- 




FIG. 4: The master diagram for the CDS of Eqs. (|ll2p . r] =0.2. There exists no CDS below the black 
solid line (here a <0, i.e. the vacuum of ([T]) is unstable). Black solid curve corresponds to q =0. Gray 
curve divides the regions, where the + and — branches of ([9]) exist. Crosses (circles) correspond to 
the +(— ) branch for ?=0.01. 



It is of interest to compare the described structure of a master diagram with the results 



presented in Refs. j2J, |25|, where the concept of "dissipative soliton resonance" has been 
raised. It is clear from the above, that this concept is interrelated with our concept of "energy 
scalable" CDS. That is a dissipative soliton resonance corresponds to an isogain asymptote of 
^ oo in Figs. Ispi (see also |9|, ll7|). For instance, the approach of 2J] based on the method 
of moments predicts the resonance condition 1/c ~ 3.16 for a =0 and x =0 versus 1/c :^1.5 in 
Fig. [31 This quantitative disagreement disappears, when the resonance condition is found on 
basis of the numerical simulations 25| instead of the approximate method of j2J| . In agreement 
with 2A1 l25|l . the asymptotic value of 1/c increases with x >0 (Fig- HI and |9(]) and decreases 
with X <0 Pi- One can conclude, that our approximate approach is more general and accurate 
than that of Ref. 2J]. Moreover, the energy scalable CDSs found on basis of our approach are 



easily traceable because they are lying on three-dimensional surfaces of reduced parametrical 
space of the CGLE. 



B. CGLE with a perfectly saturable SAM 



Such a type of SAM corresponds to a semiconductor saturable absorber, which is extensively 
used for obtaining the CDS in both solid-state and fiber oscillators j26-28|. If the CSP width 
(usually, few picoseconds) excesses the SAM relaxation time (hundreds of femtoseconds), 
the SAM function can be written in the form 



T.[\a{z,t)f] 



^nP {z, t) 

1 + KP{z,iy 



(18) 



where k = Tr/EgS is the inverse saturation power. Eg is the saturation energy fluency, 5* is the 
mode area, x =0 will be assumed below. 

The above considered method results in js], [H, 29 
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a . 2 3 



kP (0) = — = — 2(l-0-c± 

dn _ a (A^ - n^) (s^ - n^) 
m ~ 3^^^"^T]^T7A^^' 

-S^ = — + 1 - ^ + c, 
/i 3/i 

p w ^ H (A - w ) 



(19) 



where q = a/ fi, c = a^/PfiK and T = (2 — c)^ — 4<^ (2 — <^ + c). 

It should be noted, that our technique based on the adiabatic approximation and the regu- 
larization of dQ/dt j^, 0, 16 -2^, 29| is analogous to that of Ref. [30]. However, an approximate 
integration in the spectral domain allows us the further reduction of parametric space dimension 
and the construction of physically meaningful master diagrams, which make the CDS properties 
to be easily traceable. 

The CDS energy is 



E 



P (t) dt 



duj 
Zn 



(^2 _ A2 - 7//3c;) arctanh (I) 
AS 



(20) 



The normalizations of the previous subsection are useful in this case with the replacements 
K — )■ fiK, ( ^ K. Then the dimensionless CDS energy is 



r^2 



— A^ — c) arctanh (A/: 
AH 



(21) 



The master diagram following from Eq. 0211) is presented in Fig. O The CDS is two- 
parametric in this case. The structure of diagram is similar to that described in the previous 
subsection, but there are two important differences: i) there is no maximum value of 1/c- 
parameter with growth of E^ ii) — branch is not energy-scalable within all range of existence. 

Eq. fl2Tl) allows obtaining the practical important expression for the threshold energy of 
-(-branch of CDS providing its stabilization against an excitation of the vacuum of Eq. ([1]). 
Such an expression correspond to <^ =0 in Eq. (|T9|) with -|- sign before a/T (black solid curve 
in Fig. |5]): 



E. 



th 



6(2-c) 3(2 + c(c-2)) 
c Vl-c(c-3) 



arctanh 



6c (2 - c) 



2^/1 -c(c- 3) 



(22) 



Eq. ( 122|) allows estimating the threshold GDD from the desired level of the CDS energy 
and oscillator parameters. Table [T] based on Eq. ( 122|) presents the threshold values of GDD for 
a Yb:YAG thin-disk oscillator with 1 MHz repetition rate, 50-i-lOO W average output power, 
and the different semiconductor saturable absorbers (courtesy of R.Graf, Max-Planck-Institut 
fiir Quantenoptik, Garching and Dr.A.Apolonski, Ludwig-Maximilians-Universitat Miinchen, 
G arching). 
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FIG. 5: The master diagram for the CDS of Eqs. (jlllSp . There exists no CDS below the black solid 
line (here a <0, i.e. the vacuum of ([1]) is unstable). Black solid curve corresponds to ? =0. Gray 
curve divides the regions, where the + and — branches of (|19p exist. Crosses (circles) correspond to 
the +(-) branches for ?=0.1 (1) and ?=0.01 (2). 

TABLE I: Parameters (courtesy of R.Graf and Dr.A.Apolonski) of a Yb:YAG thin-disk oscillator with 
10% output coupler: a =3.6x10'^ fs^, Tr=l ps. 7 =0.22 GW~^ corresponds to an airless resonator 
with 4 transits through an active medium. Beam sizes on an active medium and a semiconductor 
saturable absorber are 2.5 mm and 0.8 mm, respectively. The spectrum is centered at «1 /um. 



Set 




2A, ps-i 






Eg, /xJ/cm^ 




1 


0.024 


0.6 


50 




90 


0.002 




0.03 


0.6 


100 




90 


0.002 


2 


0.014 


0.8 


50 




35 


0.0042 




0.018 


0.8 


100 




35 


0.0042 


3 


0.0085 


1.7 


50 




61 


0.018 




0.011 


1.7 


100 




61 


0.018 



One can see from Table [J that the CDS provides the output pulse energy levels of 50-^100 
/xJ for the /9-levels, which are achievable by the multipass chirped-mirrors delay lines. The 
minimum /3 and the maximum spectral widths are provided by the largest modulation depth /i 
(Set 3 in Table [T]). Insertion of an additional SPM source (e.g., a fused silica plate) will result 
in the (3 growth, but also it will increase the CDS spectral width and, thereby, reduce the pulse 
width after compression. 

III. CONCLUSION 

The theory of chirped dissipative solitons is exposed on the basis of the approximate in- 
tegration method of the nonlinear complex Ginzburg-Landau equation. Two main types of 
nonlinearity are considered: cubic-quintic and perfectly saturable ones, but the theory is not 
confined to a functional form of nonlinearity because only its instantaneous response is required. 
The main assumption allowing the approximate solution is dominance of nondissipative factors 
over the dissipative ones. 

The advantage of the obtained CDS solutions, which have truncated spectra with concave, 
convex and concave-convex vertexes, is their reduced parametric dimension. This makes the 
CDS to be easily traceable within an extremely broad range of parameters. 
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An application of the theory under consideration to an energy-scalable mode-locked oscillator 
is presented. It is pointed, that scaling properties of CDS allow reaching the sub-mJ energies 
of sub-picosecond pulses directly from a thin-disk solid-state oscillator operating at the MHz 
repetition rate. This promise a breakdown in further high-intensity laser experiments. 

This work was supported by the Austrian Fonds zur Forderung der wissenschaftlichen 
Forschung (FWF Project No. P20293). 
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